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The mass of the W boson is measured at LEP by fully reconstructing the W boson decays. 
The measurement techniques and systematic uncertainties are presented. The current mea- 
surement of the mass of the W boson at LEP yields 80.412 ± 0.042 GeV. 



1 Introduction 

At LEP, W bosons are produced in the reaction e + e~ — ► W + W~ with the subsequent decay of 
the W's into quark pairs, qq', or a lepton and a neutrino, tv. About 40000 W pair events are 
registered by the four experiments ALEPH, DELPHI, L3 and OPAL, corresponding to a total 
luminosity of 2.8 fb _1 . 

One of the main goals of the LEP programme is to determine the mass of the W boson, 
M\v, from the reconstructed invariant mass spectra. Involved techniques are used to obtain an 
optimal statistical precision. However, in the fully hadronic channel, systematic uncertainties 
are important, like final state interactions (FSI) between the hadronically decaying W bosons. 
These uncertainties reduce the sensitivity of this channel. The recent activities concentrate on 
increasing the weight of the hadronic events to obtain a globally more precise result on the LEP 
W mass. 



2 Extraction of the W Mass 

In the semi-leptonic and fully hadronic WW decay channels the complete four-fermion final 
states, shown in Figure can be reconstructed. The most sensitive observable for the W 
mass measurement is the invariant mass of the decaying W bosons. To improve the resolution 
on this quantity, a kinematic fitting procedure is applied. The reconstructed four- momenta of 
the final state fermions are varied within their resolution and kinematic constraints, like energy- 
momentum conservation and equal invariant masses in each event, are imposed. In fully hadronic 
events the resolution is diluted due to the different pairing combinations of the final state jets. 
In general, the most probable pairings according to the kinematics are chosen. Figure shows 
an example of the invariant mass spectrum measured in the qqqq channel. 

In ivlv final states the event kinematic can not be reconstructed completely and the leptonic 
energy spectrum and a pseudo-mass are chosen as optimal W mass estimators. 

To obtain W mass and width, three basic extraction methods are used: in the Monte Carlo 
reweighting procedure, the underlying Mw and Tw values of the Monte Carlo prediction are 
varied and compared to the measured spectra; in the convolution method the theoretically 




predicted spectra are folded with resolution functions and fitted to data; in the Breit-Wigner 
method the resonance curve is fitted to data and possible bias is calibrated with Monte Carlo 
samples. All methods exploit information from several event observables, including the invariant 
masses obtained in kinematic fits with different constraints, as well as their uncertainties. 

With the current techniques an equal statistical precision of 32 MeV and 35 MeV is reached 
in the qqiv and qqqq channel, respectively. 

3 Systematic Uncertainties 

A common source of systematic uncertainty in the qqlv and qqqq channels is the theoreti- 
cal description of photon radiation that is used in the Monte Carlo simulation. In total, the 
uncertainty on initial state radiation (ISR), final state radiation (FSR) and 0{a) electroweak 
corrections amounts to 8 MeV in all channels. The current numbers are based on reweight- 
ing of Monte Carlo events according to possible differences in the theoretical prediction of the 
photon spectrum. However, recent comparisons between two different Monte Carlo generators, 
YFSWwtDand Racoon WW^l, show larger differences. The uncertainty is therefore expected to 
be underestimated and may increase to 10 — 15 MeV for the final LEP W mass measurement. 

Another common systematic error source is the description of the hadronisation of quarks. 
The uncertainty is mainly derived from the difference between the hadronisation models Pythia, 
Herwig and Ariadne^. Especially the simulated rate of heavy baryons influences the jet masses 
and therefore the derived invariant masses. A reweighting to the measured baryon rates is 
performed by Opal. Delphi also compares mixed Lorentz boosted Z decays, that are arranged to 
reproduce W-pair decay kinematics, to the hadronisation models. Since the agreement between 
the different models and data have improved, the current systematic uncertainty of 18 — 19 MeV 
may be reduced in future. 

A similar reduction of systematic uncertainty, which is not yet included in the LEP combined 
W mass value, is originating from the LEP beam energy uncertainty. The final LEP energy 
calibration has improved H with respect to the current calibration applied. It will result in a 
much smaller uncertainty of only 10 MeV, with only a small change of the average LEP beam 
energies. The LEP experiments also performed a cross-check of the LEP energy measurement 
by reconstructing the Z boson mass in radiative fermion-pair events, as shown in Figure[21 Both 
Opal and L3 obtain results in good agreement with the LEP energy calibration^. 




Figure 2: a) Invariant Z mass spectrum reconstructed in qq7 events, b) Difference between the LEP beam energy 
and the beam energies determined by Opal in an analysis of radiative fermion-pair events. 



Z decays at high centre-of-mass energies as well as in Z peak calibration data are analysed 
to test the detector simulation. Energy and angular measurement of leptons and jets in data 
are compared to Monte Carlo simulation. Possible differences are corrected and the remaining 
uncertainty on the difference is taken as systematic uncertainty, which amounts to 14 MeV on 
the combined W mass result. 

4 Final State Interactions 

A special complication in reconstructing the invariant masses appears in the fully hadronic 
channel. Hadronic FSI may introduce cross-talk between the two decaying W bosons of one 
event. Bose-Einstein Correlations (BEC) lead to an increased production of identical bosons, 
like pions and kaons, close in phase space. Colour Reconnection (CR) changes the colour flow 
between the four quarks in the final state. By reconnecting colour strings between the previously 
colour-neutral di-quark systems from each W decay, momentum is transferred between the W's 
and the hadronisation process is altered. 

As listed in Table Q BEC and CR may change the reconstructed W mass by 35 MeV and 
90 MeV, respectively. The increase of CR effects with centre-of-mass energy is taken into ac- 
count. 

At LEP, FSI effects are measured also in other observables, which are used to constrain 
the various BEC and CR models. The BEC measurements in W pairs mainly concentrate on 
charged pion production. Assuming a spherical and gaussian shaped source emitting the pions, 
the two-particle correlation function C can be written as: 

c(Q 7r7T ) = i-\eM-R 2 QD , 

where is the square of the four-momentum difference of the two pions, — (p n ,i —p-K,2) 2 - The 
parameter A is the correlation strength and R is the inverse radius of the source. When analysing 
semi-hadronic W events, it is found that the BEC between pions coming from the same W boson 
agree very well with the correlations observed in Z decays, if Z — > bb is suppressed. 

Important for the mass measurement are, however, the correlations between pions from 
different W bosons. If there are no such correlations the two-particle density in fully hadronic 
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43 


Statistical in absence of Systematics 


32 


28 
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Table 1: Systematic uncertainties in the W mass measurement. 



events can be split into three terms: 
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The first two terms on the right-hand side of the equation are the density functions for pions com- 
ing from the same W. They can be determined in semi-hadronic events. The last term describes 
the case when one pion comes from one W boson and the second from the other W boson. This 
part can be constructed from a sample of mixed semi-hadronic event: /> w+ (1)/) W (2) = p m jx W • 
If the equation holds, i.e. in absence of BE correlations between two W's, the following ratio 
and difference of densities 



D 
Ap 



2P W + 2< ix 



,w+w- 



„ww 



2p 



w 



r, W+W- 

■^rmix 



are equal to 1 and 0, respectively. Figure 0a shows the quantity D as a function of Q as 
measured by Delphi. The Delphi data are consistent with moderate BEC between pions from 
different W's. However, the combination of all LEP experiments prefers the absence of those 
correlations^, as shown in Figure©). Monte Carlo studies show that the corresponding W mass 
shift is proportional to the BEC strength, so that the systematic uncertainty on Jl% due to 
BEC can in future be reduced to 10 — 15 MeV using the direct BEC measurement. 

A similar approach is made in the reduction of the CR uncertainty. Different models^ 
predict CR, for example the Sjostrand-Khose (SK) models 1 and 2, the Herwig model, or the 
Ariadne models 1-3. CR mainly manifests in a distortion of the W mass distribution. Effects on 
particle multiplicities are too small to be detected in data^. The only other observable sensitive 
to CR is the particle flow in the regions between the quark jets. Figure 0]a shows the angular 
distribution of particles in qqqq events, after rescaling the angle of the particle to the next jet 
in such a way that each of the four jets is positioned at an integer angular value. A ratio R is 
then calculated between the particle flow of the regions that connect two jets of one W (A+B) 
and the region that connects two jets of different W's (C+D): 



R 



d<p/dN(A + B) 
d(t>/dN{C + D) 



To quantify the difference in the two regions A+B and C+D, the R distribution is integrated. 
A variable r is then constructed as the ratio of integrals Rn obtained in data and in a Monte 
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Figure 3: a) The Delphi measurement of the D parameter shows a preference for BEC, visible as an enhancement 
of D for low Q. b) BEC results from the LEP experiments and their combination. The individual results are 
normalised such that the maximal BEC effect yields a value of 1 and the absence of BEC a value of 0. 



Carlo without CR: 

r(data) = R N (data) / R N (MC witout CR) 

Figure 0)3 shows the result of the LEP experiments, using all data. In case all LEP experiments 
are combined, the SKI model with 100% reconnection would yield a value of r(SKl — 100%) = 
0.891, while a value of r(data) = 0.969 ± 0.015 is obtained in the LEP measurement. Data 
therefore excludes strong CR effects as predicted by SKI with 5.2a. The preferred values of 
the SKI model parameter ki are in the range [0.39,2.13]^. The upper boundary of ki is used 
to estimate the systematic uncertainty on M\y due to CR. Averaged over all centre-of-mass 
energies one yields a possible mass shift of 90 MeV. This shift is larger than those observed in 
other models, like Ariadne 2 or Herwig. 

A further constraint on the CR effects can be put by the mass measurement itself. An 
indication that CR can not be strong comes from the comparison of the W mass measured in 
qqqq and qqiv events: 

AM W = Mw(qqqq) - Af w (qqft/) = +22 ± 43 MeV , 

where systematic errors due to possible FSI are removed. The observed value is compatible with 
zero. 

In addition, the W mass bias due to CR can also be reduced by modifying the jet algorithms 
that are used to reconstruct the W decay products. When reducing the cone size of the jets 
or by restricting the energy range of the objects clustered to a jet, the CR mass bias can be 
reduced, as shown in Figure EJb. Recent investigations have shown that this is the case for all 
CR models used. If one now varies the cone size or the energy cut-off, the mass shift observed 
in data and for the CR models can be compared. This result can also be combined with the 
particle flow measurement since the correlations between the measurements are small^l 

The four LEP experiments foresee to perform such an alternative mass analysis with re- 
duced sensitivity to CR effects. This will enhance the weight of the qqqq channel in the LEP 
combination, which is currently only 0.10. 



DELPHI preliminary SKI curves 




Figure 4: a) Particle flow measured by L3 in the regions inside the W's (A.B) and between W's (C,D). b) W mass 
bias due to CR for different jet reconstruction methods. Restricting the jet cone size and the energy threshold 

reduces the CR effect. 



5 Results 



Figure |S]a shows the individual results obtained by the LEP experiments. With the current 
systematic errors and all correlat ions properly included, the LEP W mass measurement split 
into qqlv and qqqq channel yields^^ 

M w (qq£u) = 80.411 ± 0.032(stat.) ± 0.030(syst.) GeV 
M w (qqqq) = 80.420 ± 0.035(stat.) ± 0.101(syst.) GeV , 

with a correlation coefficient of 0.18. The combined mass value for all channels is 

Mw(fffT) = 80.412 ± 0.029(stat.) ± 0.031(syst.) GeV , 

with a good x 2 /d.o.f of 28.2/33. Including the result derived from the cross-section measurement 
at the W-pair production threshold does not change numerically the above result: 

j\4f p = 80.412 ± 0.042 GeV . 

The method of direct reconstruction is a lso well suited to measure the width of the W boson, 
r\y- A combined fit to the LEP data yields^5 

r^f p = 2.150 ± 0.068(stat.) ± 0.060(syst.) GeV , 

and agrees well with the Standard Model (SM) prediction of Tw = 2.099 using the LEP W mass 
cited above. 

As shown in Figure EJb, the LEP W mass agrees well with the measurement at pp collid- 
ers®. The combination of these direct Mw measure ment s is also in good agreement with the 
indirect determination from the other electroweak data^l The result obtained in z^N scattering 
by NUTEV, which is derived from the measurement of the electroweak mixing angle, sin 2 9 W , 
deviates from the LEP result by 2.9 a. However, there is no systematic effect found that may 
explain this difference. 
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Figure 5: a) Measurements of Mw by the four LEP experiments, using a common systematic uncertainty for FSI. 
b) LEP result compared and combined with the measurements in pp-collisions. The direct measurement of Mw 
is in good agreement with the indirect result obtained in fits to the remaining electroweak data. 



The W mass is an important parameter in the Standard Model. The precise measurement 
of M\v probes the SM at the lev el o f its radiative corrections. A comparison of the direct and 
indirect W mass determinations is shown in Figure H3a in the plane of W mass and top 
quark mass, M t , that is measured at the Tevatron. Good agreement between the two sets of 
measurements is observed. Also shown is the SM prediction for various values of the mass of 
the Higgs boson, Mh. The measurements prefer small values of Mh- 

In a supersymmetric extension of the theory, the Minimal Supersymmetric Standard Model 
(MSSM), the preferred M\v-M t area differs slightly from the SM prediction. Additional radiative 
correction terms, involving supersymmetric particles, shift the Mw-Af t band to larger values of 
Mw- The MSSM prediction^ is shown in Figure 0b. The experimental precision of the 
electroweak measurements is not accurate enough to decide between the two models. However, 
an increased precision is expected from measurements at the Tevatron, the future LHC and 
linear colliders. The electroweak data gives also motivation to continue the search for the Higgs 
boson in a low mass range, as the last missing piece of the SM, and for supersymmetric particles 
at the Tevatron and at LHC experiments. 



6 Conclusion 

The mass of the W boson is measured at LEP to M^ EP = 80.412 ± 0.042 GeV and the width to 
= 2.150 ± 0.091 GeV. It is an important contribution to the tests of the Standard Model 
and its supersymmetric extensions. Exploiting the recent results on systematic uncertainties, 
the Mw measurement is expected to reach a final accuracy of about 35 MeV, completing the 
many precision tests of the SM performed at LEP. 
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Figure 6: a) Contour curves for the direct and indirect measurement of Mw and M% compared to the SM prediction 
for different values of A/h. b) The direct Mw and M t measurements compared to the MSSM prediction. The 
exclusion limits from searches for new particles at LEP are taken into account in the calculation. The MSSM 
and SM areas overlap in the region where the mass of the SM Higgs boson is in the MSSM range, i.e. for 

M H ~ 130 GeV. 
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